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Abstract 

We construct Poisson structures for Ermakov systems, using the Er¬ 
makov invariant as the Hamiltonian. Two classes of Poisson structures 
are obtained, one of them degenerate, in which case we derive the 
Casimir functions. In some situations, the existence of Casimir func¬ 
tions can give rise to superintegrable Ermakov systems. Finally, we 
characterize the cases where linearization of the equations of motion 
is possible. 

PACS numbers: 45.20.Jj, 05.45.-a, ll.30.-j 


1 Introduction 

Ermakov systems ]T|-|| have attracted attention in the last decades due to 
their important physical applications and mathematical properties as well. 
A central mathematical property of Ermakov systems is the existence of a 
constant of motion, the Ermakov invariant. The Ermakov invariant allows 
to construct a nonlinear superposition law linking the solutions of the equa¬ 
tions of motion composing the Ermakov system Q. Ermakov systems have 
recently been of interest in diverse scenarios, such as accelerator physics ||, 
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dielectric planar waveguides [jSj. cosmological models BB. analysis of su¬ 
persymmetric families of Newtonian free damping modes 0, study of open 
fermionic systems JlCj, analysis of the propagation of electromagnetic waves 
in one-dimensional inhomogeneous media [[□]] , algebraic approach to integra- 
bility of nonlinear systems |H|, coupled linear oscillators |13|, the semiclas- 

in 


sical limit of quantum mechanics 


supersymmetric quantum mechanics 


15 


16- 18 


computation of geometrical angles and phases for nonlinear systems 
search for Noether |T9], ^(J and Lie |21], |22| symmetries, the pos- 

extension of the Ermakov system 


sible linearization of the system 
concept |^T 


23, 24 


25|-|27|, the search for additional constants of motion [f28l and 


some discretizations of Ermakov systems [£f 

The existence of a Hamiltonian or Lagrangian formulation is a central 
question for any dynamical system. Cervero and Lejarreta [[H|] have identified 
a Hamiltonian subclass of Ermakov systems and have used this Hamiltonian 
formulation as the starting point for the quantization of these systems. Lat¬ 
ter on, Haas and Goedert have extended the class of Hamiltonian Ermakov 
systems by inclusion of frequency functions depending not only on time, but 
on dynamical variables as well |[82]| . Both Hamiltonian formulations for Er¬ 


makov systems are canonical formulations, for which the Poisson bracket 
is defined in the conventional way. On the other hand, non canonical, or 
generalized Hamiltonian, or Poisson descriptions, have proven to be relevant 
in such diverse fields like magnetohydrodynamics, kinetic models in plasma 
physics, biological models, optics, quantum chromodynamics and so on |33 


There is a wide range of possibilities open when a Poisson formulation is 
available, such as nonlinear stability analysis through the energy-Casimir 
method, perturbation methods, integrability results and bifurcation proper¬ 


ties |33| 


Accordingly, in recent years there has been interest in constructing 
¥Tj, mainly for the special case of three-dimensional 


Poisson structures 
dynamical systems. 

A finite-dimensional dynamical system is said to be a generalized Hamil¬ 
tonian system when it can be cast in the form 


x^ = J^ v d v H , (jl — 1,..., N , (1) 

where sum convention is assumed and d^ = d/dx 11 . Here, H = H(x ) is the 
Hamiltonian function and J^ v = J^ v (x) is the Poisson matrix for the system. 
The Poisson matrix must be skew-symmetric, J^ v = —J v>x . Moreover, it must 
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satisfy the following set of partial differential equations, 

J pv d v J pa + J pv d v J ap + r u d v J pp = 0. 

These equations ensure that the generalized Poisson bracket, defined as 


( 2 ) 


{A,B} = d fi AJ^d v B, (3) 

for any functions A = A(x) and B = B(x), satisfies the Jacobi identity 


{A {B, C}} + { B , {C, A}} + {C, {A, B}} = 0 . 


(4) 


In fact, equations @ are a necessary and sufficient condition for the bracket 
(3;) to satisfy the Jacobi identity. Hereafter we refer to (|2j) as the “Jacobi 
identities” too. The above generalized Poisson brackets are endowed with 
all properties of conventional Poisson brackets, with the advantage of being 
applicable to more general systems. 

From the definition, we can identify the basic building blocks of any Pois¬ 
son formulation as being the Hamiltonian function and the Poisson matrix. 
If a time-independent constant of motion is known, the idea is trying to 
use it as the Hamiltonian function for the system. In other words, we can 
look (HD in a reverse way, as a set of equations for some of the components 
of the Poisson matrix. Then, for given x p and H , system (IT]) is an unde¬ 
termined linear system for the matrix elements J pv . Besides being skew 
symmetric, the Poisson matrix must comply with the Jacobi identities (2[), 
that constitute an overdetermined system of partial differential equations for 
the remaining components J pv not fixed by ©• This “deductive schema 
for constructing Poisson structures was developed in detail in the reference 
4l|, and applied to three-dimensional dynamical systems like Lotka-Volterra 


systems for three interacting populations, the Rabinovich system and the 
Rikitake dynamo model [jfG[]. The basic proposition of the present work is 
to put forward the approach of ref. [[Il[], to find new classes of Ermakov 
systems for which a Hamiltonian formalism is possible. We use the only 
time-independent constant of motion always available for Ermakov systems, 
namely the Ermakov invariant, as the Hamiltonian function and check the 
consequences of this assumption. This idea is partly inspired by the results 


of reference ||2|, where it was shown that (n + 1)—dimensional extensions of 
Ermakov systems, when restricted to the unit sphere S n , can sometimes be 
viewed as canonical Hamiltonian systems with the Ermakov invariant play¬ 
ing the role of Hamiltonian. Here, however, there is no restriction to any 
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particular submanifold, and we focus on non-canonical descriptions. Finally, 
notice that Ermakov systems are non autonomous four-dimensional dynam¬ 
ical systems, in contradistinction to the earlier works |4(| [4]J focused on 
three-dimensional models. 

The paper is organized as follows. In Section 2, we set the Hamiltonian 
for Ermakov systems as being the Ermakov invariant, and seek a Poisson 
matrix complying with the Jacobi identities. In this way, we arrive at two 
classes of Ermakov systems admitting Poisson structures, analyzed in detail 
in Section 3. In Section 4, we examine the possibility of applying linearization 
transforms to the resulting Ermakov systems with Poisson character. Section 
5 is dedicated to our final remarks. 


2 Poisson structures 

Classical Ermakov systems are commonly written in the form 


X + u 2 x = 

-\f(.y/x ), 

( 5 ) 


yx z 

y + uj 2 y = 

-^g(%/y) , 

xy z 

( 6 ) 


where / and g are arbitrary functions of the indicated variables, and u is 
an arbitrary frequency function. In most applications, u is restricted to be 
dependent on time only, but here this constraint is relaxed. 

For our purposes, polar coordinates r = (x 2 + y 2 ) 1 / 2 and 9 = arctan (y/x) 
are more appropriate. Ermakov systems in polar coordinates reads 

r — rO 2 + tv 2 r = -^F(9) , (7) 

rd + 2rd = ~^G(6), ( 8 ) 

where F and G, depending only on the angle variable, are arbitrary functions 
suitably related to / and g. The main property of Ermakov systems is that 
they always possess a constant of motion, the Ermakov invariant, 

I = l -{r 2 6) 2 + j° G(X)d\. (9) 

As manifest by (H), the existence of / is not affected by any particular de¬ 
pendence on u of the dynamical variables. Also, a little thought shows that 
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for frequency functions depending on dynamical variables we can set F = 0 
in equation ([T]) without any loss of generality. However, we keep F mainly 
for ease comparison with previous results on Ermakov systems. 

We want to reformulate our Ermakov system m as a generalized Hamil¬ 
tonian system as in ([!]), with all indices running from 1 to 4 since the Ermakov 
system is a system of two second order ordinary differential equations. The 
Poisson matrix J ,w is skew symmetric and should satisfy the following system 
of partial differential equations, 


J^d^J 23 + J^d^J 31 + J^d^J 12 = 0 , 

J^d^J 24 + J^d^J 41 + J fi 4 d li J 12 = 0, 

j^<V 34 + ^ 3 < V 41 + j»%j 13 = 
J^d^J 34 + J^d^J 42 + J"%J 23 = 0 . 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


If time-independent, the Hamiltonian H is a constant of motion. The only 
time-independent constant of motion always available for Ermakov systems, 
no matter are the functions F, G and u, is the Ermakov invariant. Hence it 
is natural to define 


(14) 


H — I 


and see the consequences. 

The choice of coordinates x M is a matter of convenience. Here, we choose 


(x 1 , x 2 , x 3 , x 4 ) = (r, 9, u, v ), 


(15) 



(16) 


9 = v/r 2 , 

u = — uj 2 r + {y 2 + F(6))/r 3 , 
v = —G(6)/r 2 . 


r — u , 


(17) 

(18) 

(19) 

( 20 ) 


Using the Ermakov invariant as the Hamiltonian, 



( 21 ) 
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there results, from © and that 


u 

v/r 2 

—u; 2 r + (u 2 + F(9))/r 3 


J 12 G{9) + J 14 u, 

J 24 u, 

J 32 G(£) + J 34 u. 


( 22 ) 

(23) 

(24) 


Equation (|20| ) follows automatically from H — 0 and the skew symmetry of 
J^ (see |^T|| for details). 

Let us look more closely the system ( p2| - p4| ). Equation (^4|) can be viewed 
as the dehnition of tv, while equation (^) shows that setting J 12 = 0 elimi¬ 
nates G(9) from all considerations. 

This is a convenient choice, and still leads to a large class of examples. 
We found, after long calculations, that is very hard to impose the Jacobi 
identities when J 12 ^ 0. Thus, in what follows we set J 12 = 0, leaving G(9) 
arbitrary and allowing for more general classes of Ermakov systems. 

Summing up results from (p2f|24|) and our choice for J 12 , we obtain 


J 12 

= o, 

(25) 

J 14 

= u/v, 

(26) 

=-H 

to 

= l/^ 2 , 

(27) 

a ; 2 

= {v 2 + F(9))/r A + ( J 23 G{9) - J u v)/r. 

(28) 


Our goal now is to insert (|25| - [27|) into Jacobi identities and solve for the 
remaining components of the Poisson matrix. With the solution, we can 
know what are the allowable frequencies using (|28|) . 

The second Jacobi identity, equation ([□]), gives 


J 23 


u 

r 2 v 


(29) 


Inserting this and ( |2o| - [27l ) into (jT0|), we obtain 


u 


dJ 13 dJ 13 13 

+ v—^— = J 13 


u 


du 


dv 


(30) 


with solution 

Jl3 = (u) +U ^ (u ,r,0, O ' ( 31 ) 

Here, ip is an arbitrary function of the indicated arguments. Notice that we 
have included time-dependence for extra generality. 
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Substituting the already calculated components of the Poisson matrix 
into (p~3|) gives 

dJ u 8J U , 4 2 uv , , . 

u ——— + v-z— = J 3 +- ip , (32) 


du 


dv 


whose solution is 


J 34 = 


2 uv 


$ + U V (p’ r ’ 0 ’ t ) • 


(33) 


where p is a further arbitrary function, of the indicated arguments. 

The only Jacobi identity still deserving attention is the equation (0), 
which yields the consistency condition 


z z dip v dip 2 
dr r z u dt) r 


(34) 


where the prime denotes differentiation with respect to u/v. For ip — 0, 
the consistency condition (|34D is satisfied in an immediate way leaving ip 
arbitrary. For ip ^ 0, a different class of solutions arises. 

These two possibilities and the associated Poisson structures are studied 
separately. 


3 The two classes of solution 


3.1 The case -0 = 0. 


For ip — 0, the consistency condition (|34|) imposes no constraints on the 
function ip, which remains arbitrary. From the results of the last section, we 
obtain the following Poisson matrix, 


/ 


= 


[u/v) 2 
u/ (r 2 u) 

0 


(35) 


0 0 

0 0 

— (u/v) 2 —u/(r 2 v) 

-u/v —1/r 2 

where, as said previously, p = p(u/v,r,9,t). By construction, this is a 
Poisson matrix. Moreover, it is not of Lie-Poisson, affine-linear or quadratic 


V 


-up 


u/v \ 
1/r 2 
up 

0 


type, as more usual 33 


The frequency function of the associated Ermakov system follows from 

©. 

(36) 


u 


uv 


to 2 = -(v 2 + F(9)) + —G(9) - —p(u/v , r, 9 , t) 


r 3 v 
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Using the frequency function as defined in (|3 (d|) and the definitions of u and 
v in terms of the original polar coordinates, we derive the following Ermakov 
system, 


r 

rO + 2 rO 


—tG(6) + r 2 r9p 
rW 



G(9) 




(37) 

(38) 


We see that the function F disappears from all considerations. By construc¬ 
tion, (137[ - [38D is an Ermakov system having a Poisson formulation, with the 


Hamiltonian being the Ermakov invariant and the Poisson matrix given by 
(35). This is an infinite family of Ermakov systems, containing two arbitrary 
functions, G and p. Also notice that the frequencies as given by ( |35| ) can not 
be functions of time only, necessarily having a dependence on the dynamical 
variables. 

Here, the Poisson structure is degenerate, as can be readily seen by 


det(J^) = 0. (39) 

Therefore, we can obtain Casimir functions, that is, functions which Poisson 
commute with any function defined on phase space. The defining equations 
for the Casimir functions, denoted by C, are 

J^d v C = 0 . (40) 


The existence of non constant solutions is due to the degenerate character of 
the Poisson structure. 

There are systematic methods j|3], f|4| for obtaining the Casimirs, but 
here a direct approach is sufficient. Using the Poisson matrix, we find the 
following equations for the Casimirs, 


l_dC_ 
r 2 89 


u- + v- 
8u 


udC_ 
v dr 


+ mp- 


0, 

0 . 


(41) 

(42) 


The first equation of this system shows that C do depends only on the vari¬ 
ables (■ u/v,r,9,t ). Taking into consideration this information, we transform 
the equation (|42|) into 


1 8C 8C , .dC 

->-w +a jr +a * (a ' r ' e ' t) ^ = 


0 , 


(43) 





where a = u/v. The solution for ([I3|) strongly depends on the details of 
the function ip. Notice that G(9), the extra arbitrary function defining the 
Ermakov system, does not play any role it the computation of the Casimirs. 

An illuminating way to rewrite (£h|) is found by means of the change of 
coordinates 

f—l/r , 9 — 9, a = — a . (44) 

The equation for the Casimirs becomes 




dC a . pr .dC 

a— + 1/r, 9, t)— = 0 


dr 


r “ 


da 


(45) 


For 9 interpreted as the independent variable, f as coordinate and a as 
velocity, this is Liouvillc’s equation for the invariants of the equations of 
motion 

dr _ da a . ,, , 

T0 =a ' ^ = (46) 

which are also the characteristic equations for (|45|) . Notice that here the 
physical time t is a mere parameter. 

Equations ( |46|) are equivalent to the Newton equation for one-dimensio¬ 
nal motion under the force held aip/f 2 . For functions ip yielding completely 
integrablc examples of such motions, we can find all the Casimirs for the 
Poisson structure (p5|). To show a concrete example where this is possible, 
consider the case 

r 2 dV, 

v = (47) 

for an arbitrary pseudo potential V(f,t). In this case, Newton’s equation 
that follows from (|46|) are 


d 2 r 

d(P 


dV_ 

dr 


M), 


(48) 


an autonomous potential system, since 9 does not appear explicitly. As for 
any autonomous one-dimensional potential system, there is complete inte- 
grability. The constants of motion are 


£7, = i 
2 

Sb'-M- 

(49) 

C 2 = 9- 

1 r 

dX 

(50) 

V2J 

(c 1 -v(x,t)y/ 2, 


9 



respectively the energy and the additional integration constant for the equa¬ 
tions of motion. In terms of the original coordinates of the Poisson descrip¬ 
tion, the quantities ([49|-^0|) are 


«- id 


C, = 0- 


■V(l/r,t), 
rl/r d\ 


•J2 J (C,-V( 


(51) 

(52) 


which are the Casimirs for the Poisson structure fl35|) when ip is given as in 
m- When V does not contain explicitly the time, the Casimirs becomes ad¬ 
ditional constants of motion for the Ermakov system. In this case, we obtain 
a superintegrable [[Kj] class of Ermakov systems, possessing three invariants, 
namely the Ermakov invariant and the two Casimir functions. In fact, we 


can derive superintegrable Ermakov systems in all cases when fl45|) can be 
solved in closed form for the two Casimirs of the Poisson structure, and when 
<p is time-independent. 

Finally, let us examine more closely the superintegrable Ermakov systems 
with the Casimirs C\ and C 2 given in (|5T| - |52D , in the special situation for 
which dV/dt = 0. Using C 2 as in (|52| ) and the implicit function theorem, 
we obtain locally the equation for the orbits, r = r(9, C\, C 2 ). Now, using 
the Ermakov invariant, we get the angle as a function of time through the 
quadrature 

r e r 2 (A, Ci, C 2 )d\ 


t + k -J h(X,I ) 

where k is the last integration constant and 


h(8, /1 - v / 2 [l - J" G(X)dx) 


(53) 


(54) 


Locally, (|54|) gives 9 as a function of time and four integration constants, 
namely /, C 1 , C 2 and k. 

To show an example of the procedure, consider the particular case 




for which (j48|) describes a singular oscillator. Using (j47|), there results 


r 3 6 

^ r 


(55) 


(56) 
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and, from mm we obtain the following Ermakov system, 


r = — ^G(0) - r 5 0 2 , 

r A 9 

rd + 2 r6 = —- G(9) . 


(57) 

(58) 


We left the function G(9) undetermined. From the pseudo potential (|55|) and 
the orbit equation fl52|), we get 


r 


2 


2C\ 

1 + 4CK9-C 2 ) 2 ’ 


(59) 


showing a spiral motion of a particle coming arbitrarily close to the origin. 
The time-dependence of such motion is obtained from the quadrature fl53|), 
which depends on the details of the function G(9). 


3.2 The case ijj ^ 0. 

For ijj ^ 0, the consistency condition (|34| ) has a different class of solutions, 
r u / v d\ fd'ijj 1 dip 


ip = 


^ 2 (A, r, 9, t ) ' dr 




;V'(A, r, 8 ,()) + x(r, 8, tf)ip(u/v, r, 8, t ), 


(60) 


where x is an arbitrary function of the indicated arguments. Therefore, we 
obtain the Poisson structure 


J^ = 


( 


\ 


0 0 (u/v) 2 + ui/> u/v 

0 0 u/(r 2 v ) 1/r 2 

■{u/v) 2 — —u/{r 2 v) 0 mf + 2uv^/r 

—u/v —1/r 2 —uif — 2uvi\)/r 0 


\ 


(61) 

with f specified in terms of x and i/> according to (|60|) . The Poisson structure 
is non degenerate, 


det(J-) = ^ + V>)#0. 

G v 


(62) 


The following frequency functions are derived from fl28[) , 
1 


u} 2 = — {v 2 +F{9))+^-G{9)- — (ip(u/v, r, 9, t ) + 2 -^{u/v, r, 9, t) \ . (63) 

ry> a ry> \ ry* J 
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again necessarily depending on the dynamical variables. The resulting Er¬ 
makov systems are 


f = —-^G(8)+r 2 rO +2rfhli(-^,,r,e,t)'j , 

r» + m = 


(64) 

(65) 


These Ermakov systems contains three arbitrary functions, namely G, -0 and 
X■ Since the Poisson structure is non degenerate, there are no nontrivial 
Casimirs. Therefore, we found a new class of Ermakov systems that can be 
cast in a non degenerate generalized Hamiltonian form. 


4 Linearization 


Ermakov systems with frequency functions depending only on time have 
shown to be linearizable by use of 


r = l/r 


( 66 ) 


as the dependent variable and the angle 9 as the independent one. This 
change of variables is accomplished by the relation 


6 = h{9,I)/r 2 


(67) 


where h(6,I ) is defined in In terms of f, 9 , the Ermakov systems 

transforms into an one-parameter family of second-order linear ordinary dif¬ 
ferential equations depending on the value of the Ermakov invariant [23 


whenever the frequency function does not contain dynamical variables. This 
result was used for the analysis of the stability and periodicity of some Er¬ 
makov systems arising in two-layer, shallow water wave theory [[hj. The 


linearization transform (p6| - |67|) was shown to be useful also for a class of 
Ermakov systems for which the frequency function is not a mere function of 
time |24]]. This result provides an explanation for the reasons why Kepler- 
Ermakov systems KTI, a perturbation of conventional Ermakov systems, are 


linearizable. 

In view of the usefulness of the linearization transform (|64|67|) in several 
contexts, it is interesting to check if it can be also possible for our classes 
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of Ermakov systems admitting Poisson formulations. For instance, applying 
(|66| - |67| ) to the Ermakov system (|37| - |38l) , the result is 


d 2 r 
d9 2 


1 dr 
f 2 d6 




dr 1 n ^ 


( 68 ) 


For dip/dt 7 ^ 0, (|68D becomes an integro-differential equation, a possibility 
we will not consider here. Equation (^ 8 ]) is equivalent to the equation ([46|) 
determining the Casimirs for the Poisson structure. 

The term on the right hand side of (^ 8 |) has linear character if and only 

+ B(e) f + C{8 ), (69) 


de 1 r 


de 


for functions A, B and C depending only on the angle. If the assumption 

falls in the class of linearizable 


(|69| ) is satished, the Ermakov systems 

Ermakov systems discussed in [^4[. Moreover, if fl69[) is valid the character¬ 
istic equations (ff 6 |) for the Casimirs are also linear. This does not imply, 


of course, that the Casimirs may always be found in closed form when 
holds. Similar remarks apply to the linearization of our second class of Er¬ 
makov systems admitting a Poisson formulation, treated in subsection 3.2. 


5 Conclusion 

In this paper we have proposed the Ermakov invariant as the Hamiltonian 
function and we reformulate the Ermakov system as a Poisson system. The 
main difficulty is to fold a Poisson matrix reproducing the equations of motion 
and, at the same time, being compatible with the Jacobi identities. How¬ 
ever, the task can be achieved, and two classes of Poisson structures were de¬ 
rived. One of them is degenerate, thus opening the possibility of constructing 
Casimir invariants. These Casimirs, if time-independent, are also constants 
of motion. In the cases where the Casimirs are time-independent and avail¬ 
able in closed form, we obtain superintegrable Ermakov systems. A class 
of such superintegrable Ermakov systems was explicitly shown in Section 
3.1. Another, non degenerate, class of Ermakov systems admitting general¬ 
ized Hamiltonian formulation with the Ermakov invariant playing the role of 
Hamiltonian function was also found. In this latter case, no Casimirs exist. 
Both classes of Ermakov systems are specified by two arbitrary functions. 
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All these considerations apply to frequency functions having a dependence 
on the dynamical variables. Finally, the possibility of linearization of the 
equations of motion was analyzed in Section IV. Interestingly, we found that 
the determining equations for the Casimirs are linear when the associated 
Ermakov systems are linearizable through ( |66l - |67| ). 
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